A FREE-BOUNDARY PROBLEM FOR CONCRETE CARBONATION: 
RIGOROUS JUSTIFICATION OF THE VT-LAW OF PROPAGATION 

TOYOHIKO AIKI* AND ADRIAN MUNTEANt 

Abstract. We study a one-dimensional free-boundary problem describing the penetration of 
carbonation fronts (free reaction-triggered interfaces) in concrete. A couple of decades ago, it was 
observed experimentally that the penetration depth versus time curve (say s(t) vs. t) behaves like 
s(t) = Csft for sufficiently large times t > (with C a positive constant). Consequently, many 
fitting arguments solely based on this experimental law were used to predict the large-time behavior 
of carbonation fronts in real structures, a theoretical justification of the v^-law being lacking until 
now. 

The aim of this paper is to fill this gap by justifying rigorously the experimentally guessed 
asymptotic behavior. We have previously proven the upper bound s(t) < C'Vi for some constant 
C"; now we show the optimality of the rate by proving the right nontrivial lower estimate, i.e. there 
exists C" > such that s(t) > C"\/t. Additionally, we obtain weak solutions to the free-boundary 
problem for the case when the measure of the initial domain vanishes. In this way, our mathematical 
model is now allowing for the appearance of a moving carbonation front - a scenario that until was 
hard to handle from the analysis point of view. 

Key words. Free-boundary problem, concrete carbonation, large-time behavior, \/i-law of 
propagation, appearance of a carbonation front, phase change; 
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1. Introduction. 

1.1. Background. Environmental impact on concrete parts of buildings results 
in a variety of unwanted chemical and chemically-induced mechanical changes. The 
bulk of these changes leads to damaging and destabilization of the concrete itself 
or of the reinforcement embedded in the concrete. One important destabilization 
factor is the drop in pH near the steel bars induced by carbonation of the alkaline 
constituents; see for instance [HI [T2J [25] and [5D] for technical details and [TH [3] for 
an introduction to the mathematical modeling of the situation^] The destabilization 
is caused by atmospheric carbon dioxide diffusing in the dry parts and reacting in the 
wet parts of the concrete pores. The phenomenon is considered as one of the major 
processes inducing corrosion in concrete. A particular feature of carbonation is the 
formation of macroscopic sharp reaction interfaces or thin reaction layers that progress 
into the unsaturated concrete-based materials. The deeper cause for the formation of 
these patterns is not quite clear, although the major chemical and physical reasons 
seem to be known. 

Mathematically, the proposed model is a coupled system of semi-linear partial 
differential equations posed in a single ID moving domains. The moving interface 
(front position in ID) is assumed to be triggered by a fast chemical reaction - the 
carbonation reaction. Non-linear transmission conditions of Rankine-Hugoniot type 
are imposed across the inner boundary that separates the carbonated regions from 
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the uncarbonated ones. The movement of the carbonated region is determined via a 
non-local dynamics law. 

The key objective is not only to understand the movement of a macroscopic sharp 
reaction front in concrete but rather to predict the penetration depth after a sufficient 
large time. 

A couple of decades ago, it was observed experimentally that the penetration 
depth versus time curve (say s(t) vs. t) behaves like s(t) = C\Ji for sufficiently large 
times t > (with C a positive constant). Consequently, many fitting arguments 
solely based on this experimental law were used to predict the large-time behavior of 
carbonation fronts in real structures, a theoretical justification of the vi-law being 
lacking until now. 

This is the place where our paper contributes: We want to fill this gap by justifying 
rigorously the experimentally guessed asymptotic behavior. 

1.2. Basic carbonation scenario a moving one-phase approach. We study 
a one-dimensional free boundary problem system arising in the modeling of concrete 
carbonation problem. We consider that the concrete occupies the infinite interval 
(0, oo) and that there exists a sharp interface x = s(t), t > separating the car- 
bonated from the uncarbonated zone. The whole process can be seen as a solid- 
solid phase change; see the two colors in Fig. |1.1| (left). One color points out to 
CaCO-j, (carbonated phase), while the other one indicates Ca(OH)2 (uncarbonated 
phase). The zone of interest is only one of the solid phases, namely the carbon- 
ated zone. We denote it by Q S (T) and, in mathematical terms, this is defined by 
Q S (T) := {(t,x) : < t < T, < x < s(t)} for some T > 0. Throughout this paper 
u and v denote the mass concentrations of CO2 in air and water, respectively. As 



mentioned in [2], s, u and v satisfy the following system P =P(sq, uq, vq, g, h) ( 1.1 ) 



(1.7| 



(1.1) u t - (k 1 u x ) x = f(u,v) in Q S (T), 

(1.2) v t - (k 2 v x ) x = -f(u,v) mQ a (T), 

(1.3) u(t,0) =g(t),v(t,0) = h(t) forO<i<T, 

(1.4) s'{t) = i)(u(t, s(t))) for < t < T, 

(1.5) -Kiu x (t,s(t))=ip(u(t,s(t))) + s'(t)u(t,s(t)) forO<t<T, 

(1.6) -K 2 V x {t, s(t)) = s'(t)v(t, s(t)) for < t < T, 

(1.7) s(0) = s and u(0,x) = Uo,v(0,x) — v Q for < x < s Q , 

where K4 (resp. k 2 ) is a diffusion constant of CO2 in air (resp. water), f{u,v) := 
/3(7^ — u) is an effective Henry's law, where [3 and 7 are positive constants, g and h 
are given functions corresponding to boundary conditions for u and v, respectively, 
ip(r) := a\[r] + \ p for r € R describes the rate of the carbonation reaction, where p > 1 
and a is a positive constant^] sq > is the initial position of the free boundary, while 
Mo and t>o are the initial concentrations. 

First mathematical models with free boundaries for describing the concrete car- 
bonation process have been proposed by Muntean and Bohm in [TH [TB] , where the 
first mathematical results concerning the global existence and uniqueness of weak so- 
lutions as well as the stability of the solutions with respect to data and parameters 



2 The exponent p is sometimes called order of the chemical reaction, while the parameter a is 
just a proportionality constant. Its sensitivity with respect to the model output (s, u, v) has been 
studied numerically in I17| . 
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have been investigated. Recently, we have improved their results by focussing a re- 
duced free-boundary model still able to capturing the basic features of the carbonation 
process; see [5] for the reduced model and [TJ |2| for the list of the new theorems on the 
existence and uniqueness of weak solutions to P. This model is in some sense minimal: 
It includes the transport of species (diffusion) , their averaged transfer across air- water 
interfaces (the Henry law), as well as fast reaction (with an indefinitely large chemical 
compound - "the concrete"). We have used further the advantageous structure of the 
reduced model to study the large-time behavior of the penetration depths. Basically, 
we started to wonder whether the experimentally known \fi-law 

s(t) = C'Vt for t > 0, 

where C is a positive constant, is true or not |18) . Let us comment a bit on the 
context: It was shown in [22] (pp- 193-199) that the carbonation front behaves 
like a similarity solution to a one-phase Stefan-like problem |19j . Using matched- 
asymptotics techniques, the fast-reaction limit (for large Thiele moduli) done in |15) 
for a reaction-diffusion system also led to a v^-behavior of the carbonation front 
supporting experimental results from |121 113) . e.g. On the other hand, experimental 




Fig. 1.1. (Left) Typical result of the phenolphthalein test on a partially carbonated sample 
(Courtesy of Prof. Dr. Max Setzer, University of Duisburg- Essen, Germany). The dark region 
indicates the uncarbonated part, while the brighter one points out the carbonated part. The two 
regions are separated by a sharp interface moving inwards the material. In this colorimetric test, this 
macroscopic interface corresponds to a drop in pH below 10. (Right) Computed interface positions 
vs. measured penetration depths [T7f. 

results from [30] indicate that, depending on the type of the cement, a variety of t@ 
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front behaviors with /3 ^ \ are possible. Furthermore, Souplet, Fila and collaborators 
(compare |21[ [7]) have shown that, under certain conditions, non-homogeneous Stefan- 
like problem can lead to asymptotics like s(t) ~ is. Somehow, the major question 
remains: 

What is the correct asymptotics of the carbonation front propagation? 
The main result of our preliminary investigations (based on the reduced FBP) is 
reported in 4J and supports the fact that 

s(t) — >• oo as t — > oo and s(t) < C'Vt for t > 0, 

where C is a positive constant. Moreover, in this paper we establish a result on the 
lower estimate for the free boundary s as follows: For some positive constant c 

(1.8) s(t) > cVi for t > 0. 

This estimate combined with the corresponding lower one would immediately guar- 
antee the correctness of the v^-law from a mathematical modeling point of view. In 
section [2] we derive the missing lower bound. 



Note that since, generally, fci ^> k 2 and ^ u, the system (1.1 ) — ( 1.7) cannot be 
reduced to a scalar equation, where the use of Green functions representation [5J [TO] 
would very much facilitate the obtaining of non-trivial lower bounds on concentrations, 
and hence, on the free boundary velocity. Furthermore, by using the similar method as 



the one used in the proof of ( 1.8 ), we can construct a weak solution to P satisfying sq = 
0. It is worth mentioning that Fasano and Primicerio (cf. [8j, e.g.) have investigated 
a one-phase Stefan problem when the measure of the initial domain vanishes. In 
their proof the comparison principle is used in an essential manner. However, for our 
problem P we do not have any comparison theorem for the free boundary. Our idea 
here is to develop a method to obtain improved uniform estimates for solutions and 
then use these estimates to prove the existence of weak solutions for the case sq = 0. 
This program is realized in section [3] There are neither physical nor mathematical 
reasons to believe that uniqueness of weak solutions for the case Sq = would not hold. 
However, since our fixing-domain technique is not applicable anymore, the uniqueness 
seems to be difficult to prove. 

2. Large-time behavior of the free boundary. In order to give a statement 
of our result on the large-time behavior of (weak) solutions, we consider the problem 
P posed in the cylindrical domain Q(T) :— (0,T) x (0, 1). To this end, we use the 
following change of variables: 

Let 

(2.1) u(t, y) = u(t, s{t)y) and v{t, y) = v(t, s(t)y) for {t, y) € Q(T). 

Then, it holds that 

u t - ^u vv - ^yu y = f (u , v) in Q ' (T) , 

u(0,t) = g(t),v(0,t) = h(t) for < t < T, 
s'{t)=i){u{%\)) forO<<<T, 

--^-u y {t, 1) = s'(t)u(t, 1) + s'(t) for < t < T, 
s(t) 
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— %v y {t, 1) = s'(t)v(t, 1) for < t < T, 
s(t) 

s(0) = s , u(0, y) = u (y),v(0, y) = v Q (y) for < y < 1, 

where u (y) = u (s Q y) and © (y) = u (s y) for y G [0, 1]. 

For simplicity, we introduce some notations as follows: H := L 2 (0, 1), X := {z G 
ff^O, 1) : z(0) = 0}, X* is the dual space of A, 

V(T) := J L oo (0,T;i7)ni 2 (0,T; J ff 1 (0,l)) 

and 

Vo(T) :=V(T)nL 2 (0,T;X), 

and (•,•)# an d denote the usual inner product of H and the duality pairing 

between X and X*, respectively. 

First of all, we define a weak solution of P(so, n , Vq, g, h). To do this, we use a 
similar concept of weak solution as the one introduced in [T] . 

Definition 2.1. Lets be a junction on [0,T] andu, v be functions onQ s (T) for 



< T < oo, and u and v be functions defined by (2.1). We call that a triplet {s, u, v} 

is a weak solution of P on [0, T] if the conditions (SI) ~ (S5) hold: 

(Sl)se W 1 ' oo (0,T) ri/is> on [0,T] ; (it,u) G (W 1,2 (0, T; X*)nV r (r)nL 00 (Q(T))) 2 

u-g,v-h£ L 2 (0,T;X), u(0) = w and u(0) = w - 

s'(i) = ^(it(t,s(t)) /or a.e. f G [0,T] and s(Q) = s . 



f/££) / (u t ,z) x dt+ / -^u y Zydydt+ / —(«(•, 1) + 

Jo JQ(T) s Jo s 



(/(it, u) H yu y )zdydt for z G Vb(T). 

Q(T) s 



f55J / (v t ,z) x dt+ / -^v y Zydydt+ / — l)z(-, l)di 

Jo JQ(T) s Jo s 

s' 

(-/(«, t>) H yv y )zdydt for z G Vb(T). 



Q(T) 

Moreover, let s be a function on [0, oo), and it and u oe functions on Q s :— {{t, x)\t > 
0,0 < x < s(t)}. We say i/tai {s,u,v} is a weak solution of P on [0,oo) if for any 
T > the triplet {s, u, v} is a weak solution of P on [0, T] . 

Before recalling our results concerning the global existence and uniqueness of weak 
solutions to P on the time interval [0,T], T > 0, we give the following assumptions 
for the involved data and model parameters: 

(Al) /(it, v) = (3(yv — u) for any (it, v) G R 2 where /3 and 7 are positive constants. 

(A2) g,he W/ O ' c 2 ([0,oo)) nL°°(0,oo), and g > and h > on (0,oo). 

(A3) u G L°°(0, s ) and v G L°°(0, s ) with u > and u > on (0, s ). 

Theorem 2.2. (cf. Theorems 1.1 and 1.2, Lemma 4.IJ) If (Al) ~ (A3) hold, 
then P has one and only one weak nonnegative solution on [0, 00). 
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The next theorem is the main result of this paper. 

Theorem 2.3. If g(t) — g*, h(t) = for t G [0,oo), where g* and are 
positive constants with 7/1* = g*, and (Al) and (A3) hold, then there exists a positive 
constant c such that 

s(t) > cVt for t > 0. 



The proof of Theorem 2.3 relies on three technical lemmas. We give these auxiliary 
results in the following. 

Lemma 2.4. (cf. Q Lemma 3.3]) If (Al) ~ (A3) hold, then a weak solution 
{s,u,v} on [0,oo) satisfies 

S (t) 1 nt 

Xu(t)dx + -\s(t)\ 2 + Kl / u(T,s(T))dT 



2 

,(t) ,t ,s(r) 

xv(t)dx + k 2 / / v x (r,x)dxdT 



Jo 

so 



xuodx + — I so 1 2 + / xvodx + Ki I g(r)dr fort>0. 
2 Jo Jo 

Proof. Let T > 0. In (S4) we can take z(t) = s 2 (t)y for t G [0, T] so that we have 

(u t (t),s 2 (t)y)xdt+ [ K 1 u y {t)dydt+ [ s'(t)s(t)(u(t, 1) + l)dt 
JQ(T) Jo 



(2.2) =/ ( yf (u(t),v{t)) + ^-ly 2 u y (t))s\t)dydt. 



*'(')„•.. 

Q(T) s {t) 

Here, we note that 

i-T 

(u t (t),s 2 (t)y) x dt 

= -2/ f u(t) S '(t)s(t)ydydt+ [ u{T)s 2 {T)ydy- [ u(0)s 2 (0)ydy 
Jo Jo Jo Jo 

f T / >s(t) s'(t) f s{T) f s ° 

(2.3) ——2 / / u(t) xdxdt + / u{T)xdx — I u$xdx; 

Jo Jo S W Jo Jo 

(2.4) f KiUy(t)dydt = m [ (u(t,l) - g(t))dt; 
Jq(t) Jo 

(2.5) [ T S -^(u(t, 1) + l)s 2 (t)dt = C s'(t)s(t)u(t, l)dt + \(s 2 (T) - s 2 ); 



and 



s(t) 

r s '(t) 

-y 2 u y {t)s 2 {t)dydt 



Q(T) s (t) 

s'{t)s(t)yu(t)dydt+ [ s' (t)s(t)u(t,l)dt 
Q(T) Jo 

(2.6) =-2/ S -^-u{t)dxdt+ { s'(t)s(t)xu(t,l)dt. 

Jq(t) Jo 
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By substituting (2.3 1 < 

s(T) 



(2.6| into (2.2) we see that 



xu{T)dx + -s 2 (T) + Ki 

1 o 



xugdx 



«1 



g(t)dt 



u(t,s(t))dt 

f(u, v)xdxdt. 

Qs(T) 



Similarly, it follows from (S5) with z(t) = s 2 (t)y that 

s(T) 



Xv(T)dx + K2 



QAT) 



v x (t, x)dxdt 



xvodx 



Qs(T) 



f(u, v)xdxdt. 



Adding these two equations leads to the end of the proof of this lemma. □ 

Before starting off to providing a proof for the main result of the paper (Theorem 



2.3), we wish to point out in Lemma 2.5 and in Lemma 2.6 below that our free- 



boundary problem allows for positive and uniformly bounded concentrations, and 
also, that an energy-like inequality holds. 

Lemma 2.5. (cf. E] Lemma 3.2]) Assume (Al) ~ (A3) hold, take positive 
numbers g* and h* satisfying uq < g* , v < h* on [0, So], g < g* , h < h* on [0,oo) 
and g* = 7/1*, and let {s,u, v} be a weak solution of P on [0, 00). Then it holds that 



0<u<g*,0<v<h* onQ s 



Lemma 2.6. (cf. 01 Lemma 3.4]) Under the same assumptions as in Theorem 



2.3 a weak solution {s,u, v} of P on [0,oo) satisfies 

-J \ u {t)-g if \ 2 dx+ 1 -J^ \ v (t)-h*\ 2 d X +-J | S '(r)| 1+2 /^T 

rt r s ( T ) r f Mr) 

+K-L / / \u x (r)\ 2 dxdT + 7K2 / / \v x (r)\ 2 dxdT 
Jo Jo Jo Jo 

1 f s ° 7 f s ° 

<o / K - g*\ 2 dx + - / \v -K\ 2 dx 



2 Jo 1 u a " 2„ 
+ / s'(r)^\g,\ 2 +g, + ^\K\ 2 )dT fort>0. 



2.1. Proof of Theorem |2.3[ . In this section, we give the proof of our main 
result. 

Proof. Let {s,u, v} be a weak solution of P on [0,oo), and g* and h* be positive 
constants defined in Lemma 2.5 First, Lemma 2.6 implies that 

,t r s(r) 



< 



< 



7«2 
1 

2 
1 



Jo 



l^z (f) I dxdr 
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uo - g*\ 2 + 7ko - K\ 2 )dx + / s{T){-\gA 2 +g* + ],\h*\ 2 )dT 



\ua - g*\ 2 + l\ v o - K\ 2 )dx 



+ ( M 2 +9* + 7,\h*\ 2 ) (*(*) ~ so) for * > 0. 
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Hence, there is a positive constant depending on uo, vq, <?*, /i* and so such that 

ns(r) 
lUa^pttedT < Ci + Cis(t) for i > 0. 



Next, on account of Lemma 2.4 we see that 

,(t) 1 rt 

xu(t)dx + -\s(t)\ 2 + K\ I u(t,s(t))cIt 
2 Jo 

s(t) /■« /.s(r) 

xv(t)dx + K2 / / v x (T,x)dxdr 
Jo Jo 

xugdx + — I so | 2 + / xvodx + n\ I g(r)dr 
2 Jo Jo 

> K%g*t for t > 0. 

Here, we note that 



i(i,s(*)) = ( — ) 1/p fort>0. 
a 



Then, by putting M — max{g* , h*} we obtain 

r s(t) 



K X g*t < 2M J xdx + -\s(t)\ 2 + (s'(r)) 1/p dr 

+^ 2 ( / H 2 cfedT) 1/2 ( / dxdr) 1/2 for * > 0. 



It is clear that 



y v K " - a i/p 



a 

(2.7) < ^Ki5*i + C 2 s(t) fori>0, 

where C2 is some positive constant, and 

k 2 ( / \v x \ 2 dxd T ) 1/2 { f dxdT) 1 ' 2 

<K 2 ci /2 (i + s (t))^ 2 ty 2 s(t)^ 2 

< -jKi9*y + C 3 (s(t) + s(t) 2 ) fori>0, 

where C3 is some positive constant. 

From the above inequalities we can get 

~ki«M < (M + - + C 3 )| s (t)| 2 + (C 2 + C 3 )a(t) for i > 0. 

Now, let i > 1. In this case we see that 

-Ki5,t < (M+ - + C 3 )|s(i)| 2 + C 4 |s(£)| 2 + i« lfl r»t for t > 0, 
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where C4 is some positive constant. Thus it holds that 
In case < t < 1, we have sov^ < s(t). 

( k \ 1/2 
Therefore, by putting v$ = min{so, ( 4(m+i+c 3 +c 4 ) ) } we conc lude that 

u Vt < s(t) for t > 0. 

□ 

3. Appearance of a moving carbonation front The case s = 0. The 

aim of this section is to prove a result concerning the existence of weak solutions to 
P for the case sq = 0. This is the case when the free boundary starts off moving 
precisely from the outer boundary [exposed to CO2]. Before giving the statement of 
the theorem, we denote for simplicity 

C ((0,T];X) = {ze C([0,T] :X):z = 0on [OA) for some 5 Z > 0}. 

Theorem 3.1. Let T > 0, and g and h be functions on [0,T] satisfying g,h G 
V^ 1,2 (0,T) and g(t) > So > and h > for t G [0, T], where go is a given positive 
constant. Then under (Al) there exists a triplet {s,u,v} of functions such that s £ 
W^°°(0,T), s(0) = 0, s(t) > for t G (0,T], u,v G L°°{Q(T)), u - g,v - h G 
L 2 (0,T;X), u,tJGC7((0 j r];fr), fi, u G W^((0, T]; X*), 

(3.1) s'(t) = 1)) /or a.e. te [0,T], 



/ (u t ,z) x dt+ / -^u y Zydydt+ / — z(-,l)dt 
Jo JQ(T) s Jo s 



-^UyZydydt + I 
lQ(T) " 

(3.2) =/ (/(u,«) - -yu y )zdydt for z G C ((0,T];X), 

JQ(T) S 



T f K2 

(v u z) x dt+ / -^VyZydydt 
JQ(t) s 

(3.3) =-/ (/(fi,«) + -|^ v )zdydt /orze C ((0,T];X), 

JQ(T) S 



where u and v are functions defined by (2.1) 



Proof. First, let {son} be a sequence satisfying son > for each n and so n — > as 
n — > 00 and put uq u = 5(0) and «on = h(0) on [0, s 0ll ]. Then, Theorem |2 . 2 1 guarantees 
that P(so n ,uo n ,vo n , g,h) has a unique weak solution {s n ,u n ,v n } on [0,T]. Here, we 
denote by u n and v n the functions defined by ( |2.1[ ) with s = s„, u = u n and t; = v„ 
for each n. Since we can take positive constants g* and h* such that g < g* and 



/i < h* on [0,T], uon < 3* and Uo« < h* on [0, son] for n and g* = 7/1*, Lemma 2.5 
implies that 

(3.4) < u n < g*, < v n < h* on Q S „(T) for any n. 
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By (S3) and this shows that |s^(t)| < ip(g*) for t 6 [0, T] and n so that the set 
{s n } is bounded in W 1 ' 00 (0,T). Clearly, there exists a positive constant L\ such that 
< s n (t) < Li for t € [0, T] and n. 

Next, the following estimate is a direct consequence of [TJ Lemma 4.2]: For each 

n 

«i / / |u na; | 2 dxdr + n 2 / / |w n2: | 2 rfa;(iT 
Jo Jo Jo Jo 

< 2(C) + 1) f f " ( Vn(r) - ff(r)| 2 + |u„(r) - h{ T )\ 2 )dxd T 



10 Jo 

- -1 I 1 s n (r)(\f(g(r) 7 h(r))\ 2 + \ 9 r(r)\ 2 + \h T {r)\ 2 )dr 



+ / s' n (T)(3\g(r)\ 2 + \g(r)\ + \h(r)\ 2 )dr for t € [0,T], 



where Cf := /3j. Because of the boundedness of {s n } and (3.4) there exists a positive 
constant M 2 such that 

,(t) /-T /.a n (r) 

\u nx \ 2 dxdr + / / |w„ x | 2 da;dT < M2 for n. 



Then, easily, we can obtain that {u ny } and {w ny } are bounded in L 2 {Q(T)). 

From now on we provide the estimate from below for the free boundary as follows. 
To do so from Lemma [23] it follows that 



.(*) 1 ft 

Kigot < I xu n (t)dx + -\s n {t)\ 2 + kx I u„(t, s n (r))dT 

xv n (t)dx + k 2 / / v nx (r,x)dxdT 
Jo Jo 

=■ J\n{t) + J2n(t) + J 3 «(t) + Ji n (t) + Jhn{t) for i > and n. 

Here, it is obvious that 

Jm(t) + J in (t) <\{g* + h*)\ Sn (t)\ 2 \ 

. > ,1/2, ,,,,„ > lor t € 0,T and n. 

J 5n (t) < K2 M^ /2 (ts n (t))^ 2 J 1 J 



Similarly to (2.7 1, by using ( |l.4[ ) we observe that 

J 3 n(t) < f K{Tt ,P dT < ^'/T S ^) 1/P fOT 1 e t°' T ] and n - 

From the above inequalities we have 
nig t 

< 9 * +h * + l s n {t) 2 + KlTl v T S ^) lh + ^(M 2 T)^ 2 s n (ty/ 2 

2 ar' p 

< ^+h* + l L2r + + K2(M2T) l/ 2L l/ 2 ^ Sn{tr 

-. M z s n {ty for t E [0,T] and n, 
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where 

fi := min{l/p, 1/2} 

so that 

(3.5) s n (t) > vxt 1 ^ for t G [0,7] and n, 

where v\ is a positive constant independent of n. 

As next step, we wish to estimate the time derivative of u n . Let 8 > and 
T) G L 2 (<5,T; A). Then (S4) implies that 

T 

(u nt (t),r)(t))xdt\ 

s 

^ I I* J £rAu y {t),r ]y {t)) H dt\ + \ [ T (<Mu n (t,l) + ^-Mt,l)dt\ 



5 



+ | / (f(u n (t),v n (t)),r)(t)) H dt\ + \ [ S ^{yu ny {t),i 1 {t)) H dt\ 

J 8 J 8 s nW 

= : An + An + An + -An- 



Obviously, on account of (3.5| it holds that 



hn < hi ^ t2/fi \u ny (t) \ H \Vy{t)\Hdt 

- '^2^/Jl\ nn y\L 2 (S,T-H)\'ny{ t )\L 2 (S,T-M)i 

An < "^Jl/ji (\Un\L*(0,T;Hi-(p,l)) + T 1/2 ) \rj\ L 2 (S,T;X) \ 

hn < P^h* + g^T 1 / 2 ^^^; 

T , i>(g*) , , , , 

hn < - £^ \Uny\L2(0,T-,H)\V\L2(5,T;X) toi n. 

Hence, the set {u n t} and {v n t} are bounded in L 2 (5,T; X*) for each 5 > 0. 

From these estimates we can take a subsequence {rij} C {n} satisfying s n , — > s 
weakly* in W 1 ' 00 ^, T) and C([0, T}), and u rl] -> u and w nj S weakly* in L°°(Q(T)) 
and weakly in L 2 (0, T; JT^O, 1)), in C([5, T];iJ) and weakly in W^ 2 (S, T; X*) for each 
5 > as j oo, where s € M /1,oo (0, T), w,v G L°°(Q(T)), u-g,v-h£ L 2 {0,T;X) 
and 

u, € G C((0, T]; if ) n Wf D ' c 2 ((0, T];X*). 



By (3.51 we have s(t) > for t > 0. Also, s(0) = 0. 



In order to complete the proof of the Theorem, it is necessary to show that (3.2) 



(3.3) and (3.1| hold. Let r) £ C ((0,T];X). Then rj = on [0,<5] for some S > 0. By 



taking z — rj in (S4) we infer that 
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+ 



(u njt (t),ri(t))xdt 



(f(u n At),V nj (t)) 



' ')v(t,l)dt 

yu„ ]y {t),r](t)) H dt for j. 



Sn 3 (t) 

s nj (t) 



Elementary calculations yield: 

(3.6) u nj {-,l) u(;l) in L 4 (5,T) as j -)• oo 

so that 



t „> 



n 3 W 



Sn, (*) 

Moreover, we can obtain s' r 



S '(t) 
S(t) 



u(t, l)r](t, l)dt as j 



s' in i 4 (5, T) as j — > oo and 



(-^jyu v (t),r](t)) H dt as j 



Therefore, we can prove that (3.2) holds. Similarly, (3.3 1 is valid. Finally, by (3.6 1 
we get (3.1 1. Thus the proof of this theorem has been finished. □ 

A practical comment. Before using in the engineering practice the \/t infor- 
mation for forecasting purposes, the practitioner should be aware of the fact that its 
validity is closely related to the validity of the underlying free-boundary model P. 
Relying on our working experience with such FBPs for carbonation (based on |17j . 
e.g.), we can say that P captures well accelerated carbonation tests, but it may not be 
suitable for predicting the evolution of carbonation scenarios under natural exposure 
conditions. 
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